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Abstract 

We demonstrate and study the enhancement of intermediate-field two-photon absorption by 
shaped femtosecond pulses having spectral phases of antisymmetric nature. The intermediate-field 
regime corresponds to pulse intensities, where the two-photon absorption is coherently induced by 
the weak-field nonresonant two-photon transitions as well as by additional resonance-mediated four- 
photon transitions. It is a regime of significant excitation yields, exceeding the weak-field yields by 
two orders of magnitudes, reaching about 10-20% population transfer. The considered antisymmet- 
ric nature is with respect to one-half of the (initial-to-final) two-photon transition frequency. The 
corresponding pulse spectrum is detuned from this frequency (the detuning direction is according 
to the system). We study in detail the coherent interference mechanism leading to the observed 
enhancement using forth-order frequency-domain perturbative analysis. We also show that, even 
though the maximal enhancement is achieved with phase patterns of perfect antisymmetry, at high 
enough intermediate-field intensities absorption enhancement beyond the transform-limited level 
is still achievable even with patterns having some degree of deviation from perfect antisymmetry. 
The degree of tolerance to deviations from perfect antisymmetry increases as the pulse intensity 
increases. The theoretical and experimental model system of the study is atomic sodium. These 
findings are of particular importance for coherent control scenarios that simultaneously involve 
multiple excitation channels. 
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271 ] . The broad spectrum 



of the femtosecond pulses photo-induces a coherent manifold of state-to-state initial-to-final 
pathways that by manipulating their interference state-to-state transition probabilities are 
controlled. Inducing constructive or destructive interferences among these pathways leads 
to either an enhancement or attenuation, respectively, of the transition probability. It is 
achieved by tailoring the shape of the femtosecond pulse. Experimentally, it is implemented 
using pulse shaping techniques that are applied in the frequency domain to control the 



spectral characteristics of the pulse 



28 



Rational design of the pulse shape is based on identifying the state-to-state interfering 
pathways and their interfering mechanisms. As the pulse shaping is conducted in the fre- 
quency domain, rational pulse design is most powerful when the potoexcitation picture is 
available in the frequency domain. This is possible once the excitation is validly described by 
the time-dependent perturbation theory of finite order, following a proper Fourier transfor- 
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2lj, where the multiphoton process is described by 



the time-dependent perturbation theory of the lowest nonvanishing order. For N-photon 
process it is the Nth perturbative order. The simplicity of this theoretical description is 
the attractiveness of the weak-field regime. Still, it has a considerable downside of being 
limited to low excitation yields, typically up to only ~0.1% population transfer. Recently, 
we have extended theoretically and experimentally the rational pulse design approach to 
the intermediate-field regime 26|, |27J], where the excitation yields exceed by two orders of 
magnitude the weak-field yields, reaching 10-20% population transfer. In this regime, the 
perturbative description of the multiphoton process needs to additionally include also the 
next non- vanishing order beyond the lowest one, i.e. orders N and N+2 for N-photon process. 
Specifically, we have studied the process of atomic two-photon absorption [2i 



2j. In 



the weak-field regime it is described by second-order perturbation theory, while in the 
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intermediate-field regime it is described by forth-order perturbation theory considering both 
the second and forth perturbative orders. Physically, the corresponding weak-field absorp- 
tion involves initial-to-final nonresonant excitation pathways that are composed of two ab- 
sorbed photons, while the intermediate-field absorption also involves additional resonance- 
mediated pathways that are composed of three absorbed photons and one emitted photon. 
Consequently as already has been shown a decade ago §j], the weak-field two-photo ab- 
sorption is maximized by the transform-limited (TL) having zero (relative) phase across the 
whole spectrum as well as by any shaped pulse having a spectral phase that is antisymmet- 
ric around one-half of the two-photon transition frequency Uf g . On the other hand, as we 
have recently shown in short 26|], in the intermediate-field regime the shaped pulses with 
antisymmetric phase patterns actually enhance the absorption beyond the level induced by 
the TL pulse. It occurs when the center of the pulse spectrum is detuned from Uf g /2, either 
to the red or to the blue according to the system. 

The enhancement of the intermediate-field two-photon absorption by shaped pulses with 
spectral phase patterns of antisymmetric nature is the subject of the present theoretical 
and experimental work. The coherent interference mechanism leading to the observed en- 
hancement is studied and analyzed in detail using forth-order frequency-domain perturbative 
analysis. We also show that, even though the maximal enhancement is achieved with phase 
patterns of perfect antisymmetry, at high-enough intermediate-field intensities absorption 
enhancement beyond the TL level is still achievable even with patterns having some degree 
of deviation from perfect antisymmetry. The corresponding degree of tolerance to deviations 
from perfect antisymmetry increases as the pulse intensity increases. The theoretical and 
experimental model system of the study is atomic sodium (Na). Beyond the basic scien- 
tific understanding, these findings are also of significance for coherent control scenarios that 
simultaneously involve multiple excitation channels. For example, the tolerance to devi- 
ation from perfect antisymmetry allows to keep the excitation of the two-photon channel 
sufficiently high, while adjusting the excitation of the other channels. 

The paper is organized as follows. Section II gives first the frequency-domain theoret- 
ical description of the intermediate-field femtosecond two-photon absorption process and 
then refers to the case when the absorption is induced by shaped pulses having phase pat- 
terns of antisymmetric nature. Section III presents the corresponding experimental and 
theoretical-numerical results for the intermediate-field two-photon absorption in Na, and 
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Sec. IV analyzes and discusses them in detail. The paper ends with conclusions in Sec. V. 



II. THEORETICAL DESCRIPTION 

A. Intermediate-field description of femtosecond two-photon absorption 

The atomic femtosecond two-photon absorption under consideration is from a ground 
state \g) to an excited state |/), which are coupled via a manifold of states \n) having the 
proper symmetry. The pulse spectrum is such that all the \g)-\n) and \f)-\n) couplings are 
non-resonant, except for the \ f)-\n r ) resonant coupling. In other words, the spectral ampli- 
tude is zero at all the \g)-\n) and |/)-|n) transition frequencies, i.e., \E{u gn )\ = \E(uf n )\ = 0, 
except for the \ f)-\n r ) transition frequency, i.e., \E{ujf nr )\ ^ 0. The corresponding excitation 
scheme is shown in Fig. [H 

Within the present intermediate-field regime, the time-dependent (complex) amplitude 
Af of state |/), following irradiation with a (shaped) temporal electric field e(t), can be 
validly described by 4' h -order time-dependent perturbation theory. Generally, it includes 
non-vanishing contributions from both the 2 nd and 4 th perturbative orders: 

A f = Af + Af. (1) 

The perturbative description allows a transformation into a frequency- domain, where the 
spectral field of the pulse E(u) = \E(u) \ exp [i$(u;)] is given as the Fourier transform of e(t), 
with |i£(a>)| and $(o;) being, respectively, the spectral amplitude and phase of frequency u>. 
For the unshaped transform-limited (TL) pulse, 3>(a>) = for any uj. We also define the 
normalized spectral field E(u) = E(uj)/ \Eq\ = E(u) exp [z$(o;)] that represents the pulse 
shape, where \E Q \ is the peak spectral amplitude. This allows to clearly distinguish in the 
expressions given below between the dependence on the pulse intensity and the dependence 
on the pulse shape. The maximal spectral intensity Jo is proportional to \Eq\ 2 (Iq oc |-Eo| 2 ) 
and correspond to different temporal peak intensities Itl of the transform-limited (TL) 
pulse. 

As shown before for the weak-field regime 8], the 2 nd -order amplitude A^ is given by 



A 



(2) 

/ ih 2 



l-lEof A^\u f9 ) , (2) 



/oo 
E(uj)E(VL — uj)duj , (3) 
-oo 



where Uf g is the |^)-|/) transition frequency and fij g is the corresponding real effective non- 
resonant two-photon coupling. Equations (j2J) and ([3]) reflect the fact that A\ coherently 
interferes all the non-resonant two-photon pathways from \g) to |/) of any combination of 
two absorbed photons with frequencies u) and uj' = ujf g — oo, i.e., having their frequency sum 
equal to Uf g . Several such two-photon pathways are shown schematically in Fig. [TJ 
The 4^-order amplitude term is given by 



4 (4) - 1 \F I 4 



a (on— res) . a [near— res) 



(4) 

A (on-res) = ^ A & {u f g ) A^\<S) , (5) 

/oo i 
d8-A^{uj fg -5)A^ R \8) , (6) 
-oo " 

where A^ 2 \Q,) is defined in Eq. and 

A {R \AQ) = A {non - resR) {AQ) + A {resR \AQ) , (7) 

/oo 
E(oj + AVL)E*(uo)duj , (8) 
-oo 

A iresR)^ m) = |^J 2 J zvrj g(^ /nr+ Afi) J E*( W/nr ) 

/■ +0 ° 1 ~ ~ 1 

-pi dd'-E(uj fnr + An-5')E*(uj fnr -5') . (9) 

J — oo " 

This set of equations reflects the fact that Ap interferes all the four-photon pathways from 
\g) to \ f) of any combination of three absorbed photons and one emitted photon. Several 
typical four-photon pathways are shown schematically in Fig. [TJ 

As seen in Fig. [TJ and reflected in the above equations, each four-photon pathway is 
either on-resonant or near-resonant with either \g) or |/) with a corresponding detuning of 
5. Accordingly, A^ [Eq.(jT])] has two contributions: A^ n res ^ [Eq.flSJ)] coherently interferes 
all the on-resonant pathways having 5 = 0, while Af* res ^ [Eq. (l5]) ] coherently interferes 
all the near- resonant pathways having 5^0. The on-resonant pathways are excluded 
from the integration in A^ ear ~ res ^ by the Cauchy's principle value operator p. Each of the 
amplitudes A^p~ res ^ and A ( j iear ~ res " > is divided into a sequence of two two-photon parts: (i) 
a part that interferes all the non-resonant two-photon transitions with a frequency sum of 
Q = ujf g — 5 and (ii) a part that interferes all the Raman transitions with a frequency 
difference of AQ = 5. The border line between these two parts is detuned by 5 from |/) 
or \g) according to whether, respectively, the two-photon part [(i)] or the Raman part [(ii)] 
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comes first. These two parts are expressed, respectively, by the parameterized amplitudes 
A^{Q) [Eq.©] andA( R )(Afi) [Eq.©]. 

The Raman transitions are either of non-resonant nature or of resonance-mediated nature 
via \n r ). The non-resonant Raman transitions are interfered in the parameterized amplitude 
^inon-resR) [EqJS] , with ^ and jjjjr being the \g)-\g) and |/)-|/) real effective non-resonant 
Raman couplings due to all the non-resonantly coupled states \n). The resonance-mediated 
Raman transitions are interfered in the parameterized amplitude A^ resR ^ [Eqj9], with \if nr 
being the |/)-|ji r ) dipole-moment matrix element. The corresponding on-resonant (5' = 0) 
and near-resonant {5' 7^ 0) pathways are separately interfered there, with 5' being the 
detuning from \n r ) (see Fig. [T]). 

For a given physical system and a given pulse shape E(u), a non-zero Ay is proportional 
to Iq while a non-zero A^p is proportional to 1$, i.e., their ratio is proportional to Iq. 
For a given Jo, the relative magnitude and relative sign between Aj and A^ are generally 
determined by the pulse shape E(oS) and by the magnitudes and signs of the different Raman 

couplings (fJ>gg, (i 2 ff, and |/f/n r | 2 )- In the present work, for a set of intensities Io, the final 

2 (2) (4) 2 

\f) population Pf = \Af\ = Af + A\ (i.e., the degree of the two photon absorption) is 
controlled via the pulse shape E(uj). 



B. Two-photon absorption by shaped pulses with antisymmetric phase patterns 

The weak-field two-photon absorption is associated with the amplitude A^ 2 \ujf g ) [Eqs.© 
and (|3J)]. The phase associated with each \g)-\f) two-photon pathway is ^pathway = $(^) + 
$(cj/ g — uj). So, with the TL pulse all these pathways acquire zero phase ^pathway = 
and thus interfere one with the other in a fully constructive way. For a given spectrum 



\E((j) I, this leads to the maximal value of 



A f 



and to the maximal degree of weak-field 

(2) 

A^ is also 



non- resonant two-photon absorption [8] . Additionally [8J , this maximal value of 
induced by any shaped pulse having a spectral phase pattern <J>(a>) that is antisymmetric 
around one-half of the two-photon transition frequency Uf g /2, i.e., satisfying the relation 
$(u;) = — &(ujf g — u), which also yields zero phase ^pathway = for all the two-photon 
pathways. 

We consider the phase patterns $(u;) to be divided into 2n bins, each having a spec- 
tral width of Abi n . The inversion point of the antisymmetric pattern is located between 
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the bins 2n/2 and 2n/2+l and corresponds exactly to Uf g /2. So, each antisymmetric pat- 
tern contains n bin pairs such that Q m (ujf g /2 + £) = ~~ ^m{^fg/2 ~ £) for any £ satisfying 
{(m — l)Afe iri < £ < mA(, in } with m = 1 . . . n. If a bin does not have the same phase value 
as its pair partner, the phase of the corresponding pathway is non-zero (^pathway 7^ 0) and, 
hence, this pathway does not interfere constructively with the other pathways. This leads to 
the reduction of the two-photon absorption amplitude |A^ 2 ^(c<jj s )| below its maximal value. 
The larger the number of such non-matched bins is, the smaller is the degree of two-photon 
absorption. 

In the intermediate-field regime the final amplitude is contributed by both the second- 
and fourth-order terms. In the fourth-order term, the value of A^(Q = Uf g ) determines 
the contribution of the pathways that are on resonant with \g) or |/) [Eq.(J5D]. The actual 
shape of A( 2 \u}f g ± 5) is also important: following the integration in Eq.([6]) it determines 
the contribution of the pathways that are near resonant with \g) or |/). In this work the 
spectrum of the pulse is chosen such that the central spectral frequency is of red detuning 
from Uf g /2, with E{ojf g /2) ~ 0.5. Such a spectrum generally corresponds to the typical 
case, where is negligible relative to A^ in the weak-field limit and becomes comparable 
to A^ in the upper intermediate-field limit. 

We have previously shown [26| that, with a pulse spectrum that is detuned to the red 
from one-half of the two-photon transition frequency, the intermediate-field femtosecond 
two-photon absorption in atomic sodium is enhanced beyond the TL level by shaped pulses 
having the antisymmetric phase patterns. As in the weak-field regime, such patterns max- 
imize the magnitude of the second-order amplitude [A^^ = ujf g )\ to the level of the TL 
pulse. On the other hand, they lead to a fourth-order term with a much smaller magnitude 

( 4 ) (fi = uJf g )\ as compared to the TL pulse. However, the inter-term interference between 
the second- and the fourth-order amplitudes is actually of destructive nature. Hence, the 
antisymmetric phase patterns suppress this destructive interference and lead to an enhanced 
absorption as compared to the TL absorption. In the work below we demonstrate that the 
requirement for perfect antisymmetry can be relaxed to some degree. Even if the phase 
pattern contains several bin pairs that their bin members do not have the same phase mag- 
nitude and opposite sign, at sufficiently high intensities (still within the intermediate-field 
regime) such shaped pulses still enhance the two-photon absorption beyond the TL level. As 
explained below, this effect results from the destructive interference between the second- and 
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fourth-order terms, which is stronger for the TL pulse, in combination with fact that higher 
intensity leads to stronger destructiveness as the ratio between the fourth- and second-order 
terms increases with intensity (see above). Also plays a role here the interference between 
the on- and near-resonant components of the fourth-order amplitude, as they are affected 
differently by the deviation from perfect antisymmetry of the phase pattern. 



III. RESULTS 

Our model system is atomic sodium 29] with the excitation diagram shown on Fig. [TJ 
The ground state \g) corresponds to the 3s state, the excited state |/) corresponds to the 4s 
state, and the |n)-states correspond to the various p-states, with \n T ) being the 7p state. The 
two-photon transition frequency W/ )S = 00^^= 25740 cm^ 1 corresponds to two photons of 
777 nm and the transition frequency oj nr j = co , 7 P)4s =12800 cm -1 corresponds to one photon 
of 781.2 nm. The corresponding inversion point of the antisymmetric phase patterns is the 
frequency w/, 9 /2 = u;4 Si 3 S /2=12870 cm -1 . 

Experimentally, vapor of atomic sodium is produced in a static vacuum chamber at 
350°C with 10-Torr Ar buffer gas. It is irradiated at a 1-kHz repetition rate with amplified 
phase-shaped linearly-polarized femtosecond Gaussian laser pulses. The central spectral 
wavelengths is Ao=780 nm, the corresponding bandwidth (FWHM) is 5 nm, and the TL pulse 
duration is ~180 fs. The laser pulses undergo shaping in a 4/ optical setup incorporating 
a pixelated liquid- crystal spatial light phase modulator in the Fourier plane of the shaping 
setup |28(. The experiment is conducted with different pulse energies. Upon focusing, the 
corresponding temporal peak intensity of the transform-limited (TL) pulse at the peak of 
the spatial beam profile j(p™f de ~P eak ) ran g es from 5xl0 8 to 7xl0 10 W/cm 2 . Following the 
interaction with a pulse, the Na population excited to the 4s state radiatively decays to 
the lower 3p state. The fluorescence emitted in the decay of the 3p state to the 3s ground 
state serves as the relative measure for the excited 4s population Pf = P 4s . It is optically 
measured at 90° to the laser beam propagation direction using a spectrometer coupled to a 
time-gated camera system. The measured signal results from an integration over the spatial 
beam profile. 

We have divided the pulse spectral region around one-half of the two-photon transition fre- 
quency u>4 S> 3 S /2 into 2n=24 bins, with each bin having a width of Afej n =6.15cm _1 . The phase 
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values of the n=12 left-hand bins have been randomized and then inverted and assigned to 
the right-hand bins to obtain the antisymmetric pattern. Our findings are demonstrated by 
measuring the two-photon absorption induced by 2500 shaped pulses with different random 
anti-symmetric phase patterns at pulse intensities ranging from the weak-field regime to 
the upper limit of the intermediate-field regime. The results are presented as histograms 
(distributions) of the ratio between the two-photon absorption induced by the shaped pulse 
and the two-absorption induced by the TL pulse. We refer to this ratio as the enhancement 
factor. 

First, Figsf5](la) and [2](2a) present, respectively, theoretical and experimental 
enhancement-factor histograms for the set of 2500 shaped pulses having perfectly anti- 
symmetric phase patterns. The theoretical results have been calculated numerically using 
Eqs. ([I])-©, with each histogram corresponding to a different (single-valued) intensity Jo- 
Each of the spatially-integrated experimental histograms corresponds to a different pulse 
energy, corresponding to the indicated TL peak intensity at the peak of the spatial beam 
profile. The weak-field histograms are located around the enhancement value of one. As the 
intensity of the pulse increases, the histogram shifts to higher enhancement values and gets 
broader. The experimental histograms are generally broader than the theoretical ones due 
to the integration over the spatial beam profile and due to the experimental noise. The ex- 
perimental signal-to-noise ratio is directly reflected in the width of the weak-field histogram 
[Fig. E(2a), thin solid line]. 

Next, theoretical and experimental absorption results are presented for phase patterns 
that are not perfectly antisymmetric around uj^ s ^ s /2. For this purpose, in each of the random 
antisymmetric phase patterns several of the central 12 bins (out of total 24 bins) have been 
modified to be equal to their bin-pair partner, such that $(u;/ 5 /2 — £) = $(u;/ s /2 + £). 
Figs. [2](lb),(2b) and Figs. |5(lc),(2c) show the results for the almost-antisymmetric phase 
patterns with, respectively, one and two bins that have been modified. As compared to 
the perfectly-antisymmetric set, the weak-field histograms for these cases are broader and 
located below the TL value, i.e., below the enhancement factor of one. As the pulse intensity 
increases and deviates from the weak-field regime, the histogram shifts to higher values and 
an increased fraction of the pulses exceeds the TL absorption. In other words, at high 
enough intensity the absorption enhancement beyond the TL absorption is also achieved by 
phase patterns that deviate from perfect antisymmetry. 
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For the perfectly antisymmetric patterns [Fig.[2](la) and (2a)], the histograms peak for the 
maximal evaluated intensity appear at the enhancement factor of 2.2 and 1.6 for the theoreti- 
cal and experimental results, respectively. At the corresponding intensity of 2.7x 10 10 W/cm 2 , 
all the numerically evaluated phase patterns enhance the two-photon absorption above the 
factor of 1.5. For the experimental results, with a TL peak intensity of 2.9xl0 10 W/cm 2 
at the peak of the spatial beam profile, 98% of the patterns exceed the TL pulse. For the 
case of one modified bin [Fig. [2t4b) and (2b)], the weak-field histograms are peaked near 
the value of 0.9, i.e., below the TL level, for both calculated and measured results. As the 
intensity increases, 99.5% of the calculated and 93% of the measured phase patterns induce 
an absorption level above the TL level. The histograms are located, respectively, around 
the enhancement factors of 1.9 and 1.4. For the case of two modified bins [Fig. [^lc) and 
(2c)], the weak-field histograms are located around the value of 0.7. As the pulse intensity 
increases, eventually 86% of the calculated and 75% of the measured phase patterns induce 
an absorption level above the TL level. The corresponding resulted histograms are located 
around the values of 1.5 and 1.2 respectively. 

An example of a typical phase pattern leading to a high enhancement factor is presented 
in Figj3j The corresponding calculated two-photon absorption for this perfectly antisym- 
metric pattern case at the intensity of 2.7xl0 10 W/cm 2 exceeds the TL absorption by the 
factor of 3, corresponding to final 4s population of 15%. The population value is obtained 
by the exact solution of the time-dependent Srodinger equation using fourth-order Runge- 
Kutta propagation method [27] . With one and two modified bins, this phase pattern leads to 
enhancement factors of 2 and 1.35, respectively. As can be seen, there are clear phase jumps 
around the bins corresponding to the resonant transition frequencies of the excitation, u>4 St3s 
and uJip^g. Previously, step-like phase patterns have been shown to lead to the enhancement 
of resonance-mediated multiphoton absorption in the weak-field regime [ll,Q]. The corre- 
sponding mechanism is the conversion of interferences among near-resonant pathways from 
a destructive nature (for the TL pulse) into a constructive one (for the shaped pulse with 
the phase step). 
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IV. DISCUSSION 



In the discussion below we analyze the effect of the phase patterns that enhance the 
two-photon absorption in the intermediate-field regime using the perturbative description of 
Eq.(pT])-([2]). We analyze first the functional dependence of the parameterized amplitudes 
and A( R > on the detuning S from the two-photon transition frequency tOf g , which reflects the 
interference between the two-photon pathways of the corresponding type, and as it is shown 
below has a crucial effect on the two-photon absorption. Then we discuss the product 
of these amplitudes as it contributes to the fourth-order perturbative term. The aim of 
this analysis is to clarify the interference mechanism leading to the enhancement of the two- 
photon absorption for the perfectly antisymmetric and modified phase patterns as compared 
to the TL pulse. 

The amplitudes and A^ R > are shown on Fig|4]as a functions of the normalized detuning 
5/ Alj, with Au being the pulse spectral bandwidth. The real and imaginary parts of the 
complex amplitudes are shown separately on the different panels: panels FigJH 1(a) and 
1(b) for the real and imaginary part of A^ respectively. The Raman amplitude A^ is 
shown on FigJ3]2(a), 2(b), and 3(a), 3(b). In order to simplify the analysis we discuss first 
only a non-resonant part of A^ R \ i.e. A^ = A^ non - resR ^ [FigfJ 2(a) and 2(b)]. Next, we 
consider both the non-resonant and resonant parts: A^ R ' = j\i non - resR ) _|_ A( resR ) [FigH|3(a) 
and 3(b)]. Due to the p-integration with 1/5 weighting factor in the near-resonant term 
j^jiear res) .t^j], the mostly dominating region is the region with small values of 5/Auj. 
It is marked schematically on the corresponding panels. 

The term A( 2 \Q) coherently integrates all the non-resonant two-photon transition am- 
plitudes E{ijj)E{uj is ^ s — 5 — to) contributed by all the possible pairs of photons with a fre- 
quency sum of Q = uJ4 Sy 3 S — 5. The TL pulse induces fully constructive interferences among 

(2) 

all these two-photon pathways for any Q, the corresponding value of A T ' L {Q, = u^^s — 5) 
is actually the maximal one (positive and real) for any given 5 as it is drawn on FigJH 
(la), (lb). As can be seen from Eq.©, is a self-convolution of the corresponding 



spectrum 



. Thus, with a Gaussian spectrum around lo , it is peaked at £l pea k = 2co>o> 
i.e., at Speak = ^4s,3s — 2co>o. For loq < cl>4 Si3s /2 (Ao =780 nm) we obtain 5 pea k > and 
A^ (co>4 S; 3 S — 5) monotonically increases around 5 = upon the negative-to-positive increase 
of 5, i.e., A (2 \uj 4s , 3s - \5\) < AW(uo AS)3s + \8\) for small \8\; 
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The non-resonant Raman term A( n ° n - res R) coherently integrates all the non-resonant 
Raman amplitudes E(u+8)E*(u) contributed by all the possible photon pairs with frequency 
difference of 8. It is generally symmetric around 8 = 0, and obtains at that point its maximal 



E(u) 



For the case of the TL pulse, 



value which depends only on the spectral amplitude 
the amplitude contributed by any pair of photons is real and positive and so is the resulting 
total amplitude Aj^ 1 resR ^ [see Figj4](2a),(2b)]. The resonant Raman term A^^ 1 ^ integrates 
all the on-resonant Raman transition amplitudes in its first term as well as near-resonant 
Raman transition amplitudes in the second term. The amplitude corresponding to the 
photon pathways of the first type originates from the pairs of photons that are detuned 
by 8 from the |/)-state but are on- resonant with |n r )-state (8' = 0) for any possible 8: 
E(ujf nr + 8) E* (u f nr ) . Following the Cauchy principle value theorem, for the TL pulse the 
corresponding amplitude contributes to the imaginary part of A^ R ' as given in Eq.([9]). The 
second term integrates the amplitudes contributed by the photon pairs that are detuned by 
8' 7^ from the |n r )-state: E(uf nr + 8 — 8')E*{ojf nr — 8'). These amplitudes are integrated 
with the weighting factor 1/8' for any possible detuning 8' and for the TL pulse contribute to 
the real part of A^ R \ Overall for the TL pulse, AjQ has smooth Gaussian shape [FigJUpanels 
(3a), (3b)] which center is determined by the spectral amplitude with respect to the three 
relevant frequencies: u>f g /2, u>f nr , and u , which in our experiment correspond respectively 
to 777nm, 781. 2nm, and 780nm. 

For a phase pattern that is perfectly antisymmetric around 0^4,5,3,5, the value of 
A ai [ Usym {Vt = 00^^) is equal to the corresponding TL value A y T ' L {Vt = LU4 S>3s ) due to the 

fully constructive interference among all the involved two-photon pathways. Similar to the 

(2) 

TL pulse, the phases associated with these pathways are all zero and A ar l tisym (Q = u^^s) 
has real positive value. However, as £1 deviates from tu^s^s, the real part of A^ 2 \il) grad- 
ually reduces with comparable magnitude for positive and negative deviations. As it ap- 
pears on FigH] (la), (lb), the real part of A^l tisym {uj is ^ s — 8) has a peak at 8 = 0, and 
^ {^altisym (^,3* - |$|)} ~ ^ {^antisymi^ ^ ,3s + \5\) } for small \8\. Imaginary part of 

(2) 

A a ^ tisym {Vt) near 8 = has a smooth shape, rising from negative to positive values for 
the change from negative to positive detuning and obtains a zero value at 8 = 0. 

The real part of ^antisym R ^ [ see Fig. IU (2a), (2b)] is symmetrical around 
8 = 0, i.e. ^ {A { 2us~ym SR) (^ss - \8\)j = ^{a { ™u~2 sR \u4s,zs + \S\)j and ob- 
tains its maximal value at 8 = 0, while the imaginary part is antisymmetric 
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about 5 = 0: 3 {^i^m^V^s - = -9f {^iXs/m^V^s + H)} with 

(^4s,3s) r = 0. The total Raman term corresponding to perfectly antisym- 
metric phase pattern A antisym in general is not symmetric. However, the addition of the 
A( res R) p ar t approximately conserves the symmetry of the real part of ^\ non - resR ) near 
5 = 0. The imaginary part of added to the antisymmetric % \A^~™ sR) \ 

( R) 

gives approximate symmetry of the total Raman term A a J tisym . Both real and imagi- 
nary parts of A^ tisym have a peak at 5 = and A^ tisym {ui s ^ s ) = A^i ](uas,3s) ■ For the 
small enough values of the detuning 5, which are of major importance, we can state that 

^antisymi^s^s ~ 1^1) ~ ^-antisymi^^^s + 

As the individual contributions of the two-photon amplitudes arising from different types 
of pathways are described, for the sake of simplicity of the further discussion, we define a 
product of the two-photon and Raman amplitudes AA(5) = A^ (u^^ — 5)A( R \5). The 
fourth-order amplitude [Eq.(jlj)] is rewritten and given by 



4 {4) - 1 IF I 4 



/+oo 1 
d5-AA(5) 
-oo ° 



(10) 



In the discussion below we consider the cases of non-resonant and of non-resonant and 
resonant Raman transitions independently: AA^ nR \S) corresponds to A^ R > = j[( non - resR ) 
[Fig. (la), (lb)] and AA^(6) corresponds respectively to AW = A^ non - resR ^ + A^ resR ^ 
[Fig. [5] (2a), (2b)]. Relying on the analysis above we can summarize the properties of the 
products AA(S). For the non-resonant Raman transitions only: 

AA TL (0) = AA^ an J sym (0) 
M{AAt?l ym (-\5\)} « M{AAt?i ym (+\S\)} 

^ ^^intisym(~ 1^1)} ~ ~~ ^ | ^^antisym (+ 1^ I ) } 

and for both non-resonant and resonant Raman transitions: 

AA$(0) = AA^ tisym (0) 
ft{AAW. sy J-\5\)} « ^{AA^ tisym (+\6\)} 
^ {AA^ tisym (- \S\)} » %{AAi± ym (+\5\)} 

These properties suggest that for the antisymmetric phase patterns, approximate symmetry 
of AA an t is near 5 = leads to sufficient suppression of the principle value integral in 
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Eq. (jTO!) as the integrand factor 1/8 is antisymmetric about 5 = 0. On the other hand, for 
the TL pulse AAj. L is not symmetric and the value of the corresponding integral is not as 
small as for the antisymmetric phase pattern: 



.1 
i — 

8' 



d8^AA ( Tp{5) 



> 



-OO 

d8i-AA antisyrn {8) 



(11) 



and 



.1 



d8-„AAfl{8) 



> 



-OO 

db-?AA antisym (5) 



(12) 



The amplitude of the final state |/) is determined by the interference between the second 
and the fourth order amplitudes. The nature of the interference (constructive or destructive) 
between A^j and A^j is set by the pulse spectrum and by the non-resonant Raman couplings 



sum (/i3 s3s + /4 S 4s)- We have determined that for the case of atomic sodium excited with 
red detuned pulse this interference is destructive: the fourth order amplitude A^j suppresses 
the second order amplitude A±} . Recalling that A^ oc Iq and A^ oc Iq we write the value 
of the enhancement factor corresponding to the shaped pulse and normalized by the TL one: 

2 



R 



shaped 



As, shaped 



As,TL 



shaped ~ shaped 



I ^shaped ^ 



shaped-^0 \ 



A 



(2) 
TL 



+ 4 (4) 
-+- j± tl 



|1 + K T lIq\ 



(13) 



For the antisymmetric phase patterns K shaped =K antisym =l as A\^ tisym 



the analysis above it follows that \K, 



antisym \ 



< \)Ctl\, and ICtl,^ 



antisym 



(2) 

A T ' L) while from 
< 0, resulting in 



the value of R, 



antisym 



> 1. Eq. (fl3l) reflects the mechanism of the two-photon absorption 
enhancement by the antisymmetric phase patterns. As the pulse intensity Iq increases, the 
relative weight of Aj is increased and the enhancement effect is highly pronounced. 

Consider now the antisymmetric phase patterns with several modified bins. For the 
representative phase pattern [Fig. [3], the cases of one and two flipped bins are consid- 
ered, and A^(S), A^ R \8) and AA(8) are drawn on Figs. H] and [5] respectively. The in- 
terference of the pathways corresponding to the flipped bins with the rest of the path- 

(2) 

ways has turned to be a destructive one. This fact is reflected in the value of A\J for 
these phase patterns that is reduced and does not reach the corresponding value of the 



TL pulse: A 



(2) 



< < A y 

As,2bins As,lbin 4s,antisym 



| (2) 



A^ TL . The peak shape of A^(8) has also 



been changed as compared to the antisymmetric phase pattern. On the other hand, al- 
though the peak shape of A non ~ resR > has been changed, it remains symmetrical about 
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5 = and the value of A^ non resfi )(0) remains equal to that of the TL pulse as it is ac- 
tually phase-pattern-independent. Similar observations are made concerning the term: 
451(0) = 451(0) = 42^(0) = 41(0), while the shape of A™ {8) has been changed. 
Overall, the product of the amplitudes AA(5) incorporate the properties of its both com- 
ponents: AA 2bins (0) < AA lbin (0) < AA antisym (0) = AA TL (0), while the shape of AA(S) 
remains close to symmetrical to a good approximation at the mostly contributing region of 
small values of \5\. Consequently, the contribution of the on-resonant pathways (5 = 0) to 
the fourth-order amplitude is less for the modified antisymmetric phase patterns than that 
of the perfectly antisymmetric ones, but the contribution of the near-resonant pathways 
(5 7^ 0) has not been changed sufficiently, as the symmetry of the Cauchy's principle value 
integrand is approximately conserved. 

We analyze now the enhancement factor R (Eq fl3l) considering the phase patterns with 
the flipped bins. From the second perturbative order term contribution we obtain K2Uns < 

Klbin < K-antisym = ^TL = 1, while from the fourth Order One: \JC 2 bins\ < \K-\bin\ < \^antisym\ < 

\Ktl\- These results affect the values of the total amplitude A is corresponding to the weak- 

(2) 

field limit [A 4s = A\ s '} to be A 4s 

,2bins < ^45, ibin < ^4s, antisym — ^4s,tl- The corresponding 
normalized R-values in the weak-field regime are: Rantisym — 1? Rxbin — 0.7, and R2Mns — 0.4. 
However, at the higher intensities, the intensity-dependent term JCIq plays a significant role, 
increasing the value of the enhancement factor R. The mechanism of the enhancement 
follows from the interference nature between the perturbative orders: although the values of 
4s are smaller for the modified antisymmetric phase pattern, the attenuation by 4t stays 
low, close to the values of the perfectly antisymmetric pattern, as compared to the TL pulse. 
At sufficiently high intensity the resulted total population of the final state = \A4 S \ 
exceed the population corresponding to the TL pulse. At the upper intermediate-field regime 
we obtain for the phase pattern on Fig. [3] that while its enhancement factor is R a ntisym=3, 
for the modified phase patterns -Ri& in =2, and R 2 uns— 1-35. 



V. CONCLUSIONS 



In conclusion, we have identified an extended family of phase shaped-pulses that, as op- 
posed to the weak-field regime, in the intermediate-field regime enhance the two-photon 
absorption beyond the absorption level induced by the transform-limited pulse. This corre- 
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sponding phase patterns are of antisymmetric nature about one-half the two-photon transi- 
tion frequency. The phase-patterns are demonstrated to be robust to imperfections in their 
antisymmetry: even though the largest enhancement is achieved with phase patterns of per- 
fect antisymmetry, at high enough intermediate-field intensities absorption enhancement is 
still achievable even with patterns that deviate from perfect antisymmetry. The degree of 
tolerance to deviations from perfect antisymmetry increases as the pulse intensity increases. 
These findings are of importance for the development of rational selective coherent con- 
trol strategies, for reducin g th e search space when automatic feedback control optimization 



strategies are employed 



6 



30( | . and for nonlinear spectroscopy and microscopy with shaped 



femtosecond pulses. A particularly interesting scenarios is when the excitation simulta- 
neously involves multiple excitation channels 15j. Then, the tolerance to deviation from 
perfect antisymmetry allows to keep the excitation of the two-photon channel sufficiently 
high, while adjusting the excitation of the other channels. 
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FIG. 1: Excitation scheme of femtosecond two-photon absorption in the intermediate-field regime. 
The indicated levels correspond to the Na atom (not to scale). Shown are pathway examples of 
non-resonant two-photon transitions and four-photon transitions from |/) = 4s to \g) = 3s. The 
four-photon transitions involve three absorbed photons and one emitted photon in any possible 
order, and thus can be decomposed into two parts: a non-resonant two-photon transition and a 
Raman transition. The border line between these two parts can be either on-resonance or near- 
resonance with 3s or 4s (with detuning 5). The Raman transition itself can be non-resonant due 
to the np states (except for 7p) or on/near-resonance with 7p (with detuning 5'). 
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FIG. 2: Histograms for the two-photon absorption induced in Na by 2500 pulses with different 
random spectral phase patterns that are antisymmetric about i04 S ^ s /2. Upper panels: Theoretical 
fourth-order perturbative results for different (single-valued) spectral intensities corresponding to 
different TL peak intensities Iq. Lower panels: Experimental (spatially-integrated) results for 
different pulse energies corresponding to different TL peak intensities at the spatial peak profile, 
(la), (2a) - perfectly antisymmetric phase patterns; (lb), (2b), and (lc),(2c) - antisymmetric phase 
patterns with reduced antisymmetric nature: each phase pattern has one and two modified bins 
respectively. While the two-photon absorption associated with these phase patterns is lower than 
with the TL pulse in the weak-field regime, at sufficiently high intensities it excells the TL values. 
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FIG. 3: An example of highly-enhancing phase pattern that is antisymmetric about W4 S) 3 S /2. 
The corresponding population transfer induced by the two-photon absorption in the upper limit of 
the intermediate-field regime (2.7x 10 10 W/cm 2 ) is calculated to be 15%, which corresponds to an 
enhancement factor of 3 over the TL absorption. The phase pattern is characterized by step-like 
phase changes at the resonant frequencies of the system: uji s ^ s /2 (777 nm) and W7 Pi 4 S (781.2 nm). 
See the text for a detailed description. 
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FIG. 4: Theoretical results for parametrized amplitudes A^ 2 \5), A^ non - resR \5), and A^ R \5) that 
are introduced by the frequency-domain description and used for explanation of the intermediate- 
field behavior of the shaped pulses. The results for the four pulse shapes are presented: TL pulse, 
perfectly antisymmetric phase pattern given on Fig. [3l and almost antisymmetric phase patterns 
with one and two modified bins. The amplitudes are represented as a function of normalized 
detuning 5/Auj [see text] with the most contributing region of small |<5| indicated schematically by 
the dashed area. 
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FIG. 5: Theoretical results for parametrized amplitudes AA^ non ~ resR \8) and AA^ R \5) that are 
introduced by the frequency-domain description and used for explanation of the intermediate-field 
behavior of the shaped pulses. The results for the four pulse shapes are presented: TL pulse, 
perfectly antisymmetric phase pattern given on Fig. [3l and almost antisymmetric phase patterns 
with one and two modified bins. The amplitudes are represented as a function of normalized 
detuning 5/ Aw [see text] with the most contributing region of small \S\ indicated schematically by 
the dashed area. 
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